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Abstract:  The  low-frequency  sound  speed  in  a  fluid-like  kaolinite  sedi¬ 
ment  containing  air  bubbles  was  measured  using  an  acoustic  resonator  tech¬ 
nique  and  found  to  be  1 14  m/s  with  negligible  dispersion  between  100  and 
400  Hz.  The  sediment’s  void  fraction  and  bubble  size  distribution  was  deter¬ 
mined  from  volumetric  images  obtained  from  x-ray  computed  tomography 
scans.  A  simplified  version  of  Wood’s  effective  medium  model,  which  is  de¬ 
pendent  only  upon  the  ambient  pressure,  the  void  fraction,  the  sediment’s 
bulk  mass  density,  and  the  assumption  that  all  the  bubbles  are  smaller  than 
resonance  size  at  the  highest  frequency  of  interest,  described  the  measured 
sound  speed. 
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1.  Introduction 

The  acoustic  properties  of  gassy  sediments  are  important  in  a  diverse  range  of  marine  applica¬ 
tions  including  sonar  and  mine  hunting,1  civil  and  petroleum  engineering,  and  the  interpreta¬ 
tion  of  seismic  surveys,  among  others.  A  number  of  studies  have  been  performed  regarding 
scattering  from,  and  sound  propagation  and  attenuation  in,  gas-bearing  sediments.  Examples 
are  Refs.  4-7.  Despite  these  and  other  similar  studies,  the  various  models  for  sound  propagation 
in  gassy  sediments^  1  have  gone  largely  unverified.  Wilkens  and  Richardson  investigated  the 
gas-bearing  sediments  in  Eckemforde  Bay  and  observed  resonant  bubble  acoustic  effects  con¬ 
sistent  with  their  measured  bubble  size  distribution  and  the  predictions  of  the  Anderson/ 
Hampton  model/  hereafter  referred  to  as  AH.  Unfortunately,  heterogeneous  acoustic  measure¬ 
ment  techniques,  and  possibly  heterogeneous  gas  distributions  prevented  an  unambiguous 
characterization  of  the  sediment  acoustic  properties.  Recently,  Best  et  al.  obtained  wide  band 
sound  speed  and  attenuation  measurements  on  a  gassy  sediment  and  enough  environmental 
characterization  data  to  compare  the  acoustic  measurements  to  the  predictions  of  AH.  They 
found  good  agreement  between  measured  and  modeled  acoustic  attenuation  (using  a  best  fit 
bubble  size  distribution)  but  found  poor  agreement  between  measured  and  modeled  sound 
speed.  Again,  heterogeneity  of  the  bubble  distribution,  and  deviation  from  spherical  bubbles 
were  cited. 

In  the  present  study,  we  sought  to  minimize  the  experimental  complexity  and  uncer¬ 
tainty,  and  conducted  a  laboratory  experiment  with  reconstituted  kaolinite  sediments  that  con¬ 
tained  entrained  air  bubbles.  Low-frequency  (below  1  kHz)  sound  speed  measurements  were 
obtained  with  an  acoustic  resonator  technique.  A  high-frequency  (400  kHz),  bistatic  sound 
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Fig.  1.  (Color  online)  The  acoustic  resonator  and  associated  apparatus  used  to  obtain  the  low-frequency  sound  speed 
measurements. 


speed  measurement  was  attempted  but  failed  due  to  excess  attenuation.  High-resolution  x-ray 
computed  tomography  imagery  was  used  to  determine  the  bubble  size  distribution  and  overall 
void  fraction.  The  fluid-like  kaolinite  sediment,  which  had  bubbles  of  primarily  spherical  and 
spheroidal  shape,  exhibited  a  nondispersive  sound  speed  that  was  accurately  described  by  a 
simplified  version  of  Wood’s  equation.  All  measured  sound  speeds  were  at  frequencies  well 
below  the  lowest  individual  bubble  resonance  frequency  (IBRF) 

2.  Description  of  the  apparatus  and  measurement  procedure 

The  sediment  sample  was  prepared  inside  an  acrylic  core  tube  of  circular  cross  section  with 
6. 1  cm  outer  diameter,  1  mm  wall  thickness,  and  0.463  m  length.  A  latex  rubber  membrane 
(thickness=0. 1 27  mm)  was  used  to  close  the  bottom  of  the  tube.  Dried  kaolinite  was  mixed  with 
de-ionized  fresh  water  to  a  65%  water  mass  fraction.  Air  bubbles  were  entrained  during  the 
sample  preparation.  A  measurement  of  the  high  frequency  (400  kHz)  sound  speed  was  at¬ 
tempted  using  a  pulse  transmission  technique  and  apparatus  described  in  Refs.  11,  13,  and  14, 
but  it  failed  due  to  excess  attenuation  in  the  sediment.  The  low-frequency  sound  speed  was 
determined  using  the  acoustic  resonator  shown  in  Fig.  1.  The  sample-filled  core  tube  was  trans¬ 
ferred  to  the  acoustic  apparatus  and  positioned  as  shown.  The  air-sediment  interface  at  the  top 
and  a  styrofoam  block  at  the  bottom  provided  pressure  release  acoustic  boundary  conditions  to 
a  high  degree  of  approximation,  as  described  in  Ref.  15.  Acoustic  standing  waves  were  gener¬ 
ated  with  a  3-cm-diam  aluminum  piston  that  was  attached  to  an  electromechanical  shaker  via  an 
aluminum  stinger.  Band-limited  periodic  chirps  were  produced  by  a  vector  signal  analyzer 
(VS A)  and  directed  to  the  shaker  through  a  power  amplifier.  The  source  signal  was  also  digi¬ 
tized  by  the  VSA.  The  piston  was  positioned  a  few  centimeters  below  the  upper  air-sediment 
interface.  The  acoustic  pressure  inside  the  sediment  was  received  with  a  miniature  hydrophone, 
positioned  a  few  centimeters  below  the  source.  The  received  signals  were  bandpass  filtered 
( 1 0—  1 0  kHz)  and  amplified  with  a  charge  amplifier,  and  digitized  by  the  VSA.  According  to 
manufacturer  specifications,  the  source  velocity  and  receiver  pressure  responses  were  flat  to 
within  less  than  ±1  dB  in  the  experimental  frequency  range. 

Acoustic  pressure  spectra  were  then  calculated  onboard  the  VSA  by  way  of  a  transfer 
function  between  the  received  acoustic  pressure  and  the  excitation  signal  using  20  spectral  av¬ 
erages.  This  preserved  the  pressure  magnitude  and  phase  relative  to  the  excitation  signal  and 
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Table  1.  Physical  properties  used  in  the  evaluation  of  the  elastic  waveguide  model. 


Longitudinal  sound 
speed  (m/s) 

Transverse  sound  speed 
(m/s) 

Wall  material 
density  (kg/m3) 

Inner  wall 
radius  (m) 

Outer  wall 
radius  (m) 

C/=2664 

c,=  1087 

pw=  1185 

b — 0.0295 

d=0.0305 

yielded  the  coherence  function.  The  number  of  averages  and  the  amplitude  of  the  excitation 
were  chosen  to  achieve  a  near-unity  coherence  function  at  all  frequencies  of  interest,  which 
guaranteed  a  high  signal-to-noise  ratio  and  ensured  linear  bubble  behavior.1  Peak  pressures 
inside  the  sediment  were  typically  less  than  160  dB  re  1  ^tPa. 

Immediately  following  the  acoustic  resonator  measurements,  the  sediment  sample  was 
transferred  to  a  Universal  Systems  HD-500  micro-computed  x-ray  tomography  system. 
Three-dimensional  imaging  scans  (107  (im  voxel  size)  were  obtained  over  approximately  1/3 
of  the  total  sediment  volume.  A  three-dimensional  rendering  of  a  portion  of  the  scan  is  shown  in 
Mm  1 .  The  sediment  contained  relatively  evenly  dispersed  gas  bubbles  of  primarily  spherical  or 
spheroidal  bubbles.  Image  analysis  software  was  used  to  calculate  total  gas  volume  Kgas  and 
bubble  size  distribution.  Knowledge  of  the  scanned  sediment  core  tube  volume  Vtot  was  used  to 
extrapolate  global  sediment  void  fraction  ^gas^  ^tot*  The  measured  void  fraction  was  \ 
=0.0045  and  the  void  fraction  measurement  uncertainty  was  ±0.001  due  to  the  voxel  size. 

[Mm  1  Three-dimensional  rendered  views  of  subsets  of  the  micro-CT  scans  are  shown.  The 
yellow-colored  objects  are  gas  bubbles.  Neither  the  core  tube  walls  nor  the  non-gas  phases  of  the 
sediment  are  visible  (2.9  Mb)] 

The  bulk  mass  density  psed  of  gas-free  sediment  material  was  measured  with  a  Quan- 
tachrome  Penta  pycnometer.  Five  samples  taken  from  just  below  and  five  from  25  cm  below  the 
air-sediment  interface  were  measured.  No  significant  difference  was  found  between  the  upper 
and  lower  locations.  The  mean  density  was  1581.4  kg/m3  and  the  range  of  variation  of  the  five 
samples  was  ±5.6  kg/m\  The  density  measurement  uncertainty  was  1  ppt  (given  in  pycnom¬ 
eters  specifications). 

3.  Accounting  for  the  elastic  waveguide  effect 

In  a  gas-filled  acoustic  resonator,  the  walls  are  effectively  rigid.  In  a  liquid-filled  resonator, 
there  is  significant  coupling  between  the  fill-liquid  and  the  tube.  The  result  is  a  reduced  sound 
speed  relative  to  that  observed  in  an  unconfined  environment.  This  effect  is  hereafter  referred  to 
as  the  elastic  waveguide  effect.  An  exact  analytic  model  *  for  sound  propagation  in  a  finite¬ 
thickness  elastic- walled,  fluid-filled  cylindrical  tube  was  used  to  relate  the  speeds  observed  in 
the  resonator  to  the  intrinsic  sound  speed  the  material  would  exhibit  in  an  unconfined  environ¬ 
ment.  This  procedure  (and  its  validity)  is  discussed  in  Ref.  15  and  summarized  here.  Equation 
(A  1 )  of  Ref.  1 5  is  the  dispersion  relation  for  the  resonator  waveguide.  The  intrinsic  sound  speed 
c0  of  the  liquid  that  fills  the  resonator  is  an  input  parameter  to  Eq.  (A  1 ).  The  phase  speed  cph  of 
the  plane  wave  mode  is  an  output.  The  resonator  measurements  described  here  yield  the  effec¬ 
tive  phase  speed  cefr  observed  inside  the  waveguide.  The  intrinsic  sound  speed  c0  is  then  varied 
in  Eq.  (Al)  until  the  model  output  cph  matches  the  measured  value  cefr.  The  value  of  c0  that 
achieves  the  match  is  reported  as  the  intrinsic  sound  speed  of  the  sample.  The  physical  param¬ 
eters  used  in  the  elastic  waveguide  model  [Eq.  (Al),  Ref  15]  are  given  in  Tabic  1.  The  sound 
speeds  for  the  tube  walls  were  initially  taken  from  the  literature,  but  the  similar  values  that  were 
ultimately  used  were  obtained  via  resonator  calibration  with  filtered  degassed  fresh  water. 

4.  Model  of  sound  propagation  in  fluid-like  gas-bearing  sediment 

The  sediment  in  this  experiment  was  fluid-like.  It  flowed  readily  and  its  shape  was  not  main¬ 
tained  when  unsupported.  Therefore,  Wood  s  model 1  for  the  sound  speed  in  a  bubbly  liquid  was 
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Fig.  2.  (Color  online)  The  resonator  spectrum  obtained  for  the  kaolinite  sediment  The  four  lowest  resonance 
frequencies  are  identified  with  circles,  (b)  The  effective  sound  speed  ccff  inside  the  gassy  kaolinite-filled  waveguide 
is  inferred  from  these  frequencies,  where  a  least-squares  fit  yields  ceff=99.1  m/s.  After  correcting  for  the  elastic 
waveguide  effect,  the  intrinsic  sound  speed  inferred  from  the  measurements  was  c0=114  m/s  ±2  m/s,  where  the 
range  is  due  to  finite  spectral  resolution  and  the  uncertainty  in  the  resonator  length  L.  The  Wood  limit  model  [Eq.  f  1 )] 
predicted  c0=l  10-138  ms,  where  the  range  is  due  to  the  uncertainty  of  the  void  fraction. 


used,  which  requires  that  the  excitation  frequency  be  well  below  the  resonance  frequency  of  the 
largest  bubble.  For  0.002  <  *<  0.94,  and  isothermal  conditions  (a  good  approximation  for  these 
bubble  sizes)  the  model  reduces  to19 


clfgs 


Palm  +  Psedg(r'2) 

*0  ~X)Ps'd 


0) 


where  Cifgs  is  the  low-frequency  sound  speed  in  the  gassy  sediment,  Patm-  101.3  kPa  is  the  at¬ 
mospheric  pressure,  is  the  sediment  mass  density,  g=9. 8 1  m/s2  is  the  acceleration  of  gravity, 

L  is  the  length  of  the  resonator,  and  \  is  the  void  fraction.  The  terms  in  the  numerator  account 
for  the  mean  dependence  of  the  bubble  compressibility  on  the  local  hydrostatic  pressure.  In  this 
experiment,  the  largest  bubble  radius  observed  was  4.6  mm.  The  corresponding  bubble  reso¬ 
nance  frequency  is  4.91  kHz,  determined  by  Eq.  (1 )  of  Ref.  8  using  a  conservatively  low  value 
for  sediment  shear  modulus  ( G=  1 04  Pa),  also  from  Ref.  8.  The  highest  frequency  in  these  ex¬ 
periments  was  0.8  kHz,  which  rendered  the  low-frequency  approximation  valid.  Note  that  even 
if  G  was  set  to  zero,  the  lowest  bubble  resonance  frequency  would  still  be  about  six  times  higher 
than  the  highest  acoustic  excitation  frequency.  The  simplified  form  of  Eq.  (1 )  is  used  because  it 
best  illustrates  the  physical  nature  of  sound  propagation  in  near-surface  gas-bearing  fluid-like 
sediments. 


5.  Results 

The  measured  resonator  spectrum  is  shown  in  Fig.  2(a),  where  four  resonance  frequencies  are 
identified.  The  effective  sound  speed  ceff  is  inferred  from  the  slope  of  the  curve  in  Fig.  2(b).  The 
wth  resonance  frequency  is  given  by  /„=(ceff/2L)/i  and  L  is  the  length  of  the  resonator.  After 
correction  for  the  elastic  waveguide  effect  (Sec.  3),  the  intrinsic  gassy  kaolinite  sound  speed 
was  c0=114  m/s±2  m/s.  The  measurement  uncertainty  was  due  to  the  finite  spectral  reso¬ 
lution  and  the  uncertainty  of  L  (±1 .2  mm).  For  the  length,  void  fraction  and  density  reported  in 
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Table  2.  Measured  and  predicted  sound  speeds. 


Sediment  type 

Frequency  range  of 
resonator 

measurements  (Hz) 

Low-frequency 
resonator  sound 
speed  (m/s) 

Sound  speed  predicted  by 
Eq.  (1)  (m/s) 

400  kHz  pulse 
transmission  sound 
speed  (m/s) 

Gassy  kaolinite 

100-550 

114  ±2 

110-138 

n/a  due  to  excess 
attenuation 

Degassed  filtered  water 
(calibration) 

200-1000 

1494 

n/a 

1494 

Sec.  2,  the  model  [Eq.  (1 )]  predicts  a  low-frequency  gassy  sediment  sound  speed  between  1 10 
and  138  m/s.  The  range  of  the  prediction  is  due  to  the  void  fraction  measurement  uncertainty. 
Excellent  agreement  is  observed  between  the  measurement  and  the  model.  A  summary  of  the 
measured  and  predicted  sound  speeds  appear  in  Table  2. 

6.  Conclusions 

The  two  major  conclusions  of  this  work  are:  (1)  Resonator  measurements  of  gassy  sediments 
yield  the  sub-1  BRF  sound  speed  even  where  high-frequency  pulse  transmission  measurements 
fail.  (2)  A  simplified  version  of  Wood s  equation  describes  the  sound  speed  in  shallow  gas¬ 
bearing  fluid-like  sediments  in  which  the  IBRF  of  the  largest  bubble  is  well  above  the  frequency 
of  the  acoustic  excitation,  and  the  bubbles  are  evenly  distributed  throughout  the  sample.  The  gas 
phase  dominates  the  compressibility.  The  sound  speed  is  dependent  only  on  the  void  fraction, 
the  sediment  mass  density,  and  the  local  hydrostatic  pressure.  It  does  not  depend  at  all  on  the 
material  properties  of  the  sediment  grains  and  the  fluid,  other  than  through  the  bulk  sediment 
density.  This  will  cease  to  be  true  as  the  local  hydrostatic  pressure  increases,  which  in  turn 
decreases  the  compliance  contrast  between  the  gas  phase  and  the  longitudinal  and  shear  com¬ 
pliances  of  the  gas-free  phase  of  the  sediment.  The  range  of  applicability  of  this  simple  model  is 
the  subject  of  current  work. 
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